Concatenated Control Sequences based on Optimized Dynamic Decoupling 
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Two recent developments in quantum control, concatenation and optimization of pulse intervals, 
are combined to yield a strategy to suppress unwanted couplings in quantum systems to high order. 
Longitudinal relaxation and transverse dephasing can be suppressed so that systems with a small 
splitting between their energy levels can be kept isolated from their environment. The required 
number of pulses grows exponentially with the desired order but is only the square root of the 
number needed if only concatenation is used. An approximate scheme even brings the number down 
to polynomial growth. The approach is expected to be useful for quantum information and for 
high-precision nuclear magnetic resonance. 
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If one wants to measure small shifts in local energy 
levels high-precision nuclear magnetic resonance (NMR) 
is a powerful means of investigation. In order to reduce 
the line widths which limit the achievable resolution it is 
important to isolate the system under study as much as 
possible from its environment. It is common in NMR to 
do this by intricate sequences of pulses, see e.g. Ref. [J 

Similarly, one of the big hurdles in the realization of 
quantum information (QI) is to keep information stored 
without loss due to decoherence. Besides quantum error 
correction, see e.g. 0,01; dynamic decoupling (pP_) is one 
important tool to reduce decoherence 0, [a, [a, & El- 
Here, as for NMR, a sequence of pulses is used to achieve 
an effective Hamiltonian where the quantum bit (qubit) 
is as much isolated from its environment as possible. 

So in both these fascinating areas of physics the design 
of appropriate pulse sequences is of great relevance. 

The fundamental idea is to compensate the effect of a 
Hamiltonian (or a part of it) by inverting its sign by a 
Dulse. This is known since long as spin-echo technique 
which works very nicely for static magnetic fields. 
Its sequence reads pse = fr it fr- Further improvement 
is obtained by passing to the combination of two pules in 
Pcpmg = fr 7T f%r K f T [HI EH which is commonly called 
CPMG after the first letters of its inventors. The total 
duration of the cycle is t = At. Here it stands for an 
ideal, instantaneous n pulse which takes a part of the 
Hamiltonian to its negative while f T stand for the time 
evolution due to the unchanged Hamiltonian over the in- 
terval r. Both basic cylces are commonly iterated: the 
iteration of the spin echo leads to ppdd = fr (f fr) n with 
t = (n + l)r which is often called periodic dynamic de- 
coupling (PDD). The iteration of the CPMG leads to 
PiCPMG — (/t""/t) with t — 2nr which is success- 
ful in many circumstances, see for instance Refs. IH, 14- 



It is particularly robust against soft high-energy cutoffs 

Recently, two powerful extensions of DD have been 
proposed. The first is concatenation [1, [13] (CDD). By 
iterating the concatenation of pulse sequences p n +i = 
p n X p n Z p n X p n Z with po = f T Khodjasteh and Lidar 
could show that the interaction between a qubit and its 
environment can be reduced to order t n+1 . Hereby X 
stands for a tt pulse about ax, the Pauli matrix of the 
qubit in x-direction, and Z stands for a n pulse about 
(j z • The caveat of this approach is that the number of 
necessary pulses grows exponentially with a high power 
oc 4™ with the concatenation level n, thus with the power 
in t. Hence the implementation of a CDD sequence for a 
given order can be technically very demanding. 

The second extension is the optimization of the in- 
tervals between the pulses. They need not be chosen 
equidistant and this additional degree of freedom can be 
exploited for optimization [qJ leading to the UDD se- 
quence. In a dephasing model where bosons are cou- 
pled linearly to the qubit (spin-boson model) the author 
showed that the instants Sj (j £ {1, 2, ... n}) for n pulses 
are optimum if 



Sj = tsm 2 (Trj/(2(n + 1))). 



(1) 



Based on the observation that this formula does not de- 
pend on any details of the system and on the knowledge 
that most degrees of freedom behave at high tempera- 
tures like gaussian fluctuations the author concluded this 
sequence is not specific to the spin-boson model. The 
main advantage of pudd as defined by {]]) is that one 
gains an order in t with each pulse added [£j, |l5|. So no 
exponential costs arise. 

Up to n = 5 the intervals of the UDD sequence were 
found even earlier by explict general calculation, not re- 
lying on a specific model [18|. Based on numerical evi- 
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dence and on analytical evaluation of the first orders up 
to n = 14 by large-scale recursions it was subsequently 
conjectured that |T]) applies to all pure dephasing mod- 
els 15, 0|. It is not specific to the spin-boson model 
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at all. Very recently, this well-founded hypothesis could 
be turned into a theorem by mathematical proof [20| . 
The caveat of this approach so far is that it can only 
deal with pure dephasing or longitudinal relaxation, re- 
spectively. In particular for so-called low-field systems, 
where the splitting between the two energy levels of the 
qubit is small, this is a serious drawback. In such sys- 
tems, the commonly used rotating-frame approximation 
works only poorly and so dephasing and spin flips have 
to be taken into account. 

So we are faced with two interesting sequences, CDD 
and UDD, both with strong points (CDD can eliminate 
all couplings, UDD requires only linear number of pulses) 
and weak points (CDD requires exponentially large num- 
ber of pulses, UDD eliminates only pure dephasing). In 
the present work, we show that their combination leads 
to the partial combination of their advantages. 

In the following derivation we will consider ideal, in- 
stantaneous 7r pulses for simplicity. Such pulses can be 
approximated to some extent by very short pulses (but 
see also the results in Refs . \21, l22r> and, more efficiently, 
by suitably shaped pulses ij 24j . 

To derive the advantages of a concatenated UDD se- 
quence (CUDD) we first recall that a UDD sequence of 
n pulses with 7r rotations (Z) about the z direction sup- 
presses the relaxation along z [23], i.e., very little spin 
flips occur. Formally we start from the most general de- 
coherence Hamiltonian of a single qubit 



H 



E < 

■yG{0,x,y,z} 



(2) 



where <7o is the identity in the Hilbert space of the qubit 
and the A 1 are operators of the environment only. Then 
the time evolution Pu DD of the UDD sequence 



PUDD — ft-S„ Z fs n -6„-! Z ■ ■ ■ Z fs 3 -8 2 Z hz-Si Z hx 

behaves as if an effective Hamiltonian H cS acted 



(3) 



pSdd = exp(-itif eff ) + 0(at n+1 ) (4) 

where a := max 7 (||A 7 ||) is used; any operator norm || • || 
can be employed. The effective Hamiltonian H cS only 
contains powers in a x and a y which are in total even, 
besides arbitrary powers in a z . Because an even power 
in a x or <j y equals the identity and a x a y = ia z we know 



H 



off 



76{0,z} 



(5) 



i.e., it only contains dephasing terms. The order of the 
effective bath operators obviously read 



\\A 



\Af\\ = O(m & x(\\A \\,t\\A x \\ 2 ,t\\A y \\ 2 )) (6a) 

eff I 



0(max(p z ||,ip a 



(6b) 



Note that the A^ are complicated, non-linear functions 
of t due to the control sequence. 



On this level, the idea to build another DD sequence 
based on ([5]) is very appealing because it is a pure de- 
phasing Hamiltonian. The seemingly most efficient way 
is to use another UDD sequence on top of ©, but now 
about the x- or the y-direction. Unfortunately, this is not 
possible since H eS is not independent of t. Hence, if the 
UDD sequence requires the Hamiltonian to be applied for 
a given duration this cannot be done by choosing some in- 
terval t. The dependence of H cS on t is highly non-trival 
because it is non- linear and because [H cS (t), H cS (t')] f 
for t^t'. 

The situation is better if one is not aiming at arbitrary 
intervals but at commensurate intervals {tj} which can 
be made from integer multiples of a basic interval At 
like tj — rijAt with rij G N. Then one may choose t = 
At and generate exp(—irijtH cfi ) by applying the UDD 
sequence rij times. But the intervals 5j+i — Sj are in 
general not commensurate, but see below. Only its most 
basic versions, the spin-echo for n = 1, and the CPMG 
cycle for n = 2 only use commensurate intervals. This 
is already an interesting observation since the CPMG 
sequence pcpmg suppresses dephasing up to order t 3 as 
can be seen by regarding CPMG as the n = 2 case of 
either UDD 91 or CDD, see also below. 

To assess the possible gain in concatenating CPMG 
with a UDD sequence £>udd 01 m P u l ses an d duration t\ 
we estimate that the neglected term in @ is of the order 
of (ati) m+1 This is the error committed on the UDD 
level of the sequence. On the next level we consider 



PCPMG - PUDD^PUDDPUDD^PUDD 



(7) 



so that the total duration is t — At\. Then we have 
Pcpmg = 1 + R where the deviation R is estimated 
to be of the order of ||i?|| = {a cS t/8) 3 with a cS := 
max™ (P 7 ff ||). The validity of the factor 1/8 becomes 
obvious in the generalization of the CPMG to arbitrary 
level of concatenation below. 

So the total deviation -Rcpmg-udd wu ^ ^ e approxi- 
mately 



I ptot 

I^CPMG-UDD 



|| « max((a ett £/8r, (at/4) m+L ). (8) 



As expected one has to compare expressions involving 
different powers of t. The use of m > 2 is justified if 
a is much smaller than 2a. In view of © this will be 
the case if a and a eS are dominated by the ||A 7 || with 
7 6 {x,y lZ } and if |L4 Z || < ||A X || and ||A,|| < \\A y \\. 
This observation shows that already a single concatena- 
tion step on top of the optimized UDD sequences can be 
very useful. Moreover, it is possible to pass to iterated 
CPMG sequences concatenated with basic UDD building 

blocks as in p iC pMG = (Pudd^Pudd) 20 with tl = t /( 2c )' 
where c is the number of iteration cycles. Though the 
overvall order in the duration is not changed by the iter- 
ation the influence of the unwanted couplings which are 
averaged to zero becomes smaller and smaller on increas- 
ing c. This is so because the intervals between pulses 
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become smaller and smaller, see for instance Ref. [l5j for 
data for the bosonic dephasing model. 

If the order t 3 of the CPMG is insufficient further levels 
of concatenations can be added to reach higher orders. 
To see this it suffices to consider the simplemost building 
block of concatenation 



17] 



p n+ i = p n Xp n X (9) 

If the time evolution p n is governed by H^ n > — 
J2-y£to z\ a i ® ^t™^ the Magnus expansion 25 1 straight- 



forwardly implies the recursion 



A, 



= 2r„_ 
(n+l) _ A(n) 



(10a) 
(10b) 



= i(r n /2)[A^,A^] + 0(a n T^) (10c) 

between consecutive levels of concatenation where r„ is 
the total duration of p n and a n := maXyg/g.z} 11-^-7 !!■ 

Several remarks are in order: (i) The twofold con- 
catenation according to ([9]) yields the CPMG sequence: 
P2 = PqXpqXpq because X 2 = 1. (ii) If one starts with 
bath operators Aq and A z no terms in x- or y-direction 
are generated, (iii) Note that it is sufficent to consider 
the leading orders of the Magnus expansion to derive 
(flC)]) because we are only interested in the leading devia- 
tions, (iv) The number a n of X pulses grows roughly like 
2™. More precisely, it grows like a n +i — 2a n + 2(— l) n . 
Starting from ao = this implies the exponential law 
a„ = (2/3)(2"-(-l)"). 

The essential point is that in z-direction each level of 
concatenation introduces an additional factor of r n , see 



(|10cj) . So each level of concatenation suppresses decoher- 
ence by an additonal order in the total duration of the 
pulse sequence. The pure bath operator Aq does not 
affect the qubit. This kind of property was exploited in 
Refs. H to establish CDD. Hence, if we concatenate 
n times according to ([9]) we will achieve a sizable sup- 
pression of decoherence due to high powers in the total 
duration. 

More quantitatively, instead of ([8]) we have 



I ptot I 

I^cuddI 



(a ett t2- 



\n+l 



, (atZ 



-n\m+l 



(ii) 

where t is the total duration of the entire sequence. The 
factor 2~ n (™ +1 )/ 2 in the first term results from the prod- 
uct of all the t„ obeying the recursion (|10ap . The fac- 
tor 2~" in the second term results from the number 2" 
of UDD basic building blocks; each UDD interval has 
length t2~ n . Taking © into account we can conclude 
that a < a cS for short enough times. Then we ar- 
rive at the following conservative, but simple, estimate 
II-RcuddII ~ {ctt2~ m ) m+1 ) for n = m. Hence we realize 
that the combination of concatenation and interval opti- 
mization makes it indeed possible to reduce any type of 
coupling to high order. 



In assessing the gain of the advocated CUDD over the 
known CDD approach, we have to consider the number 
of necessary pulses. In order to make the t m term vanish 
the CDD approach requires 4 m pulses HE3]. The CUDD 
(concatenated UDD) sequence on concatenation level n 
made from UDD sequences with to pulses requires a n s» 
(2/3)(2 n - (-1)™) pulses of type X and to2™ pulses of 
type Z. To suppress all t m terms we choose n = m and 
the total number of pulses is about (m + 2/3)2 m . Thus 
we have achieved a substantial reduction by passing from 
4 m to (m + 2/3)2'" if we replace CDD by CUDD. So one 
factor 2 m is reduced to the linear factor to + 2/3 due to 
the UDD building blocks. Note that this pays already 
for small to > 1. This is the main general result of the 
present work. 

Unfortunately, the suppression of the dephasing com- 
ponent along er 2 could not be reduced in the same fash- 
ion from exponential to linear by employing UDD due to 
commensurability problems, see above. If one were able 
to use a UDD scheme also for the second step the num- 
ber of pulses could be brought down to the order to 2 . In 
practice, it can be promising to check out a work-around. 
The function d(x) = sin 2 (7ra;/2) on the right hand side 
of ^ (x = j/(n + 1)) has a simple shape, see Fig. 1. 
It is odd about (1/2, 1/2), it vanishes at x = 0, and it 
has vanishing slope at x = and x = 1. These fea- 
tures are reproduced by a simple third order polynomial 
^approx = —2x 3 + 3a; 2 which approximates d(x) very well 
as shown in Fig. 1. The maximum difference is only about 



0.01 so that it is to be assumed that the use of d. 
instead of d(x) is hardly noticeable in practice. 



approx 



(x) 
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x=j/(n+l) 
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FIG. 1: (color online) Right hand side d(x) — sin 2 (7ra;/2) of 
Eq. (JTJ) compared to d ap prox = —2a; 3 + 3x 2 . 



The use of d a pprox(a;) recovers commensurability. Be- 
cause x — j/(n + 1) and the third power in d apP rox(2;) 
is the highest power occurring, all instants 5j = 
rfapprox(j/(« + 1)) are integer multiples of tj (n + l) 3 [2r3|. 
This means that we have to realize (n + l) 3 intervals of 
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UDD sequences with n Z-pulses and at the n instants Sj 
another n X-pulses are needed. Hence this type of UDD 
sequence on top of a UDD sequence (approximate UDD 2 ) 
requires n(n + l) 3 + n pulses, hence only a polynomial 
law as opposed to the CUDD which still requires an ex- 
ponential number of pulses. Due to the high power of the 
polynomial (degree 4), the approximate UDD 2 pays only 
for larger values of n if (n+1) 3 JS 2™, i.e., for n 10. Fur- 
ther exploration of the approximate UDD 2 by numerical 
means is called for because of the approximation involved 
which precludes a rigorous analytic assessment. 

In summary, we have shown that UDD can be used as a 
powerful building block in more complex pulse sequences. 
Thereby the restriction of being applicable only to the 
suppression of either pure dephasing or pure relaxation 
is abolished. Three schemes are proposed. The first is 
suited for anisotropic situations where spin flips are much 
less likely than dephasing. The UDD can be inserted 
in a CPMG cycle. The CPMG can also be iterated if 
needed. This scheme yields a suppression in t 3 only due 
to the properties of the CPMG. But in view of the robust 
properties of the CPMG cycle this may be fully sufficient 
in many circumstances. 

The second, most general and flexible, scheme con- 
sists in using the UDD sequence as starting sequence of 
concatenation to achieve a concatenated UDD (CUDD). 
This scheme makes it possible to suppress decoherence 
to arbitrary order t m+1 in the total duration t if a UDD 
sequence of m pulses is concatenated m times. The nec- 
essary number of pulse scales like m2 m which is still ex- 
ponential but considerably less demanding than 4 m as in 
CDD. 

Third, an approximate scheme of two concatenated 
UDD sequences is proposed. The appoximation is needed 
on the second level to retrieve a certain degree of com- 
mensurability. The required number of pulses scales like 
to 4 so that it becomes more advantageous than CUDD 
for a sufficiently large number of pulses (m 10). 

The proposed schemes are expected to be experimen- 
tally useful for high-precision NMR where narrow line- 
widths are a prerequisite or for long-time data storage in 
realizations of quantum information devices. 
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